Dynamics of correlations in a dilute Bose gas following an interaction quench 
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We calculate the dynamics of the one and two body correlation functions in a homogeneous Bose 
gas at zero temperature following a sudden change in the interaction strength, with and without 
an underlying lattice. We focus on conceptually simple examples that highlight the features arising 
separately from interactions and band structure. In the continuum, we show that the Bogoliubov 
spectrum leads to a diffusive evolution of density correlations for short times, and ballistic at long 
times. In the lattice we find that the correlation functions develop additional oscillations. Moreover, 
the lattice dispersion induces an additional velocity scale, and some features instead propagate 
with that velocity. Finally, we discuss the time-evolution of the contact following a quench. Our 
predictions, which can be readily tested in experiments, serve as a benchmark for considering the 
dynamics of more complicated systems. 
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I. INTRODUCTION 

New higher resolution imaging techniques are allow- 
ing cold atom experiments to probe spatial correlation 
functions [l], Recently focus has turned to the dy- 
namics of these correlations following a sudden change 
in experimental parameters such as lattice depth or in- 
teraction strength [H-Q- ^ n man y ways, experiments are 
ahead of theory, as the very paradigms for thinking of 
these highly noncquilibrium experiments arc just being 
developed. Here we use the time dependent Bogoliubov 
approximation to study quenches in the weakly interact- 
ing limit. We are thereby able to organize the phenom- 
ena, for example separating out interaction effects from 
lattice effects. Although experiments to date have been 
in a strong or intermediate coupling regime and cannot 
be directly modeled by our technique, the existing data 
appears to be consistent with these organizational prin- 
ciples. 

There arc a number of fundamental issues which can 
be probed through the time dependence of correlation 
functions. For example, experiments on lattice bosons 
have been related directly to questions of causality @, Q . 
A second interesting direction is to make connections to 
the fluctuation theorems in classical statistical mechan- 
ics Q which relate equilibrium free energy differences 
to the work done. This latter physics has been stud- 
ied in mesoscopic solid state systems by looking at the 
counting statistics of electrons tunneling through quan- 
tum dots [TTj . A third direction is using these quenches 
to help understand how isolated quantum systems ap- 
proach equilibrium, and determine the timescale for the 
development of long and short range order [r34l8| . Our 
focus is in a fourth direction, namely on how quenches 
reveal the structure of the excitation spectrum. 

We study a weakly interacting Bose gas, both in free 
space and in a lattice. We use time dependent Bogoliubov 
theory to calculate how the correlation functions evolve 
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following a sudden change in the interaction strength. 
We find that in free space the density correlations spread 
diffusively at short times, crossing over to ballistic at long 
times. The speed of propagation is simply the sound 
speed. In a lattice, new features appear which propagate 
with a speed determined by the bandwidth. The speed 
of propagation crosses over to the sound speed at long 
times. We also look at single particle correlations, for 
example finding that the time for the condensate fraction 
to reach equilibrium is simply related to both the initial 
and final chemical potential. 

In addition to this "long wavelength" physics, we find 
that the short length-scale behavior of the correlation 
functions has interesting structure. Immediately follow- 
ing the quench, the two-particle correlation function de- 
velops a divergence which scales as l/|r— r | 2 . This struc- 
ture is due to the singular nature of the two-particle rela- 
tive wave function at short distances flfjj . In equilibrium, 
this singularity has attracted a lot of attention recently 
following the work of Tan, who was able to relate the 
short distance structure of the two-particle correlation 
function in a two-component Fermi gas to the internal 
energy via a quantity called the "contact" (20l - [23| . Here 
we study the dynamics of the contact (and its general- 
izations) following a quench. 

This paper is organized as follows: In Sec II we de- 
scribe our system and derive the equations of motion 
governing the dynamics of the correlation functions in 
the presence and absence of a lattice. In Sec. Ill, we cal- 
culate the dynamics of single-particle correlation func- 
tion in the continuum after two types of quenches: a 
sudden quench from a non-interacting state, and a more 
general quench between two finite values of the interac- 
tion. In Sec. IV we describe the dynamics of the two- 
particle correlation functions both in the continuum and 
in a ID lattice. We limit our discussion to length scales 
comparable to or greater than the healing length of the 
condensate. By considering simple instructive cases, we 
illuminate the role of the excitation spectrum, and com- 
pare and contrast the lattice from the continuum. In Sec 
V., we discuss the short-range physics in the two-particle 
correlation functions, and relate our results to the con- 
tact. In Sec. VI, we summarize our results. Throughout 
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this paper, wc limit our discussion to sudden quenches at 
time t = between different initial and final states. 



II. MODEL 

We start with the Hamiltonian for Bose gas interacting 
with a contact interaction: 
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The density of condensed atoms is given by uq = 
(ak=o) 2 - We perform the standard Bogoliubov trans- 
formation wherein we write the operator a^ait) = 
Uk(t)bk + vl(t)b^_ k where bk denotes the bosonic annihi- 
lation operator for the non-condensed atoms and has no 
time dependence. The coherence factors Uk{t) and Vk{t) 
are complex numbers which satisfy |it/t(i)| 2 — |vfc(i)| 2 = 1 
for all k, at all times. We always work in the regime 
where the depletion is small, and the condensate density 
is roughly equal to the total density. 
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where g = Airah 2 jm > where a is the s-wave scattering 
length and cik is the bosonic annhilation operator. In 
free space, the dispersion e(fc) = 4p-, while in a lattice, 
we use e(fc) = — 2J cos(kd), where J is the bandwidth 
and d is the lattice spacing. Throughout, our free space 
results are in 3D, while our lattice results are in ID. These 
choices of dimensionality are somewhat arbitrary. 



Ignoring all fluctuations beyond quadratic order in the 
bk operators, and neglecting the time dependence of no, 
one obtains the usual Bogoliubov Hamiltonian, which is 
readily diagonalizcd to yield the equation of motion for 
u k and v k : 
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where Ek = y^e(k)(e(k) + 2g/no) is the Bogoliubov dis- 
persion where = h 2 k 2 /2m. Here Uk(t — 0) = 



Ek 



and v k (t = 0) = -J\ (^fa _ i 



are the static solutions. Here gi and gj denote the initial 
and final interactions strengths respectively. 



A. Expressions for Correlation functions 

Here we are interested in the dynamics of two cor- 
relation functions: the non-condensed fraction (n ex = 
Sk/o( a I a fc))' an d the density-density correlation func- 
tion gW(r - r') = E q e ' q (r_r ' (PqP-q) where p q = 
at- The former is readily measured in time- 
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of-flight experiments B , w hile the latter is probed using 
Bragg spectroscopy [4 |21[ |24| , noise correlations or 
by direct in-situ measurements [5, 6]. At zero tempera- 
ture, 



>(t) = Ek(t)f 



(3) 



~(2) 

The term §g is quartic in the UkS and VkS and arises 
from correlations between the non-condensed atoms. 
These correlations become important only at extremely 
short distances S ~ a ~ 50nm, which is roughly 20 
times smaller than the typical condensate healing length 
£ ~ /im. As we are primarily concerned with S > C, 
we ignore corrections to the dynamics arising from this 
term, except in Section V. 

The first term in Eq. 2] is roughly equal to n 2 where 
n = no+n ex is the total density. The second term denotes 
correlations between the condensed and non-condensed 
atoms. In a non-interacting gas at zero temperature, 



Uk = 1 and Vk = 0, therefore n e 
independent of 5. 



and 



(2) 
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B. Self-consistency 



9s 2 \t) =n 2 + 2n n ex (t)+n Y / e lk - S (2\vk(t)\ 2 + (4) 



where 5 



l(t)Vk(t)+Uk(t)v* k (t)) 



+ 



In a more sophisticated model, we would let no be 
time-dependent. Such an elaboration is unnecessary for 
quenches where the change in the equilibrium condensate 
fraction is small. Mathematically the resulting nonlinear 
theory is much more unwieldy. 
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III. NON-CONDENSED FRACTION 
A. Quench from g; = to j; > 

We first consider the dynamics of the non-condensed 
fraction. The conceptually simplest case is to consider 
a quench from g = to some final g > 0. From 
Eq. [5J we find that Vf-(t) — —ign sm(Ekt)/Ef., thus 
nex = EkM' = (9n Q ) 2 j: k sm{E k t) 2 /El At long 
times, the sum is dominated by small values of k, where 
the dispersion is linear E^ = ck, where c = */ gn$ / m is 
the speed of sound. With this substitution, the sum is 
readily performed to yield the expression: 



»ex(t) - -^-C" 3 (l-e- 4t/w ) 



(5) 



where C, = h/ \JfrTgno is the healing length of the conden- 
sate and tmf = fr/gno is the characteristic relaxation 
time following the quench. Therefore the population of 
non-condensed particles grows linearly at short times, 
saturating on timescales t ~ tmf- At long times, the 
number of excitations created is larger than the predicted 
zero temperature equilibrium value n eXteq = £~ 3 /37r 2 . 
Given that the Bogoliubov approximation is collisionless, 
there is no reason to expect this distribution to be ther- 
mal. 



Quench from gi > to gj > 



Next, we consider a quench from some initial interac- 
tion strength gi to a final value gj. The non-condensed 
fraction is then 



(i) = ^M0)| 2 -. g/ (. 9j - 5/ )7 



, e(fc)sin(£ fc f) s 
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The first term in this expression is the number of exci- 
tations present in the system initially. The minus sign 
in front of the second term indicates that if gi > gf, the 
number of non-condensed particles decreases over time, 
and vice- versa if ^ < gf. Alternatively, the condensate 
fraction grows if gi < gf and reduces if gf > gi- 

The number of new excitations created/destroyed is 
thus proportional to gf x (<?j — gf). If gf = 0, then this 
term is identically zero and no more excitations are cre- 
ated or lost. The initially quantum depleted supcrnuid 
remains depleted, unless interactions are present to re- 
distribute quasi-momentum. 

To obtain the relevant timescale governing the re- 
distribution of quasi-momentum, we evaluate the ex- 
pression for n ex assuming a linear dispersion. Defining 
C = h/ ^JmgiTiQ as the healing length of the initial state, 
we obtain: 

n ex (t) = ^C 3 (l - 3^/(1 - g f )t/ n [lb(4fAa) - (7) 

L {U/t 2 ) 
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FIG. 1: Quasi-momentum redistribution in 3D follow- 
ing a sudden change in interactions: Plotted is the frac- 
tion of non-condensed atoms normalized to the initial excita- 
tion fraction following a sudden quench to weaker interactions 
gj /gi = 0.6. The excitation density is expressed in terms of 
£, the healing length of the initial condensate. Note that the 
analytic formula Eq. [7] (solid curve) captures the dynamics 
very well, particularly the long time limit. The characteristic 
relaxation time is set by T2 = h/gf-n, the mean-field time of 
the final state. The dashed curve is the excitation fraction at 
the final interaction strength in equilibrium. 



where g f = gf/gt, n = ^/ ^/gjgln and r 2 = h/g f n , and 
I and L denote the Bessel function of the second kind 
and the modified Struve function respectively. For small 
arguments, Io{x) — L (x) w 1 — 2x/ir + 0(x 2 ), while for 
x>l, I (x) - L (x) w 2/irx + 0{x" 3 ). 

Unlike the case of a quench from zero interactions, 
where tmf was the only relevant timescale, the dynam- 
ics here are governed by two timescales: at short times, 
the number of excitations decreases linearly in time with 
a timescale given by the geometric mean of the initial 
and final interaction strengths ti « 1/ y/WgJno . At long 
times t > T2, the number of excitations created/destroyed 
saturates to a constant proportional to — gf). 

In Fig. [T] we plot Eq. [7] Also shown is the result of 
numerically integrating Eq. [6] assuming the full Bogoli- 
ubov dispersion, showing that Eq. [7] indeed captures all 
the relevant features seen in the dynamics. 

Again the number of excitations present at long times 
after the quench is always greater than the number of ex- 
citations in the equilibrium final state. This is consistent 
with a gain in entropy. Note that in the weakly interact- 
ing regime discussed here, the fraction of excited particles 
is typically smaller than the condensate fraction by a fac- 
tor of y/noa 3 I , so the deviation of the excited fraction 
from its equilibrium value may be somewhat difficult to 
measure experimentally. The clearest (albeit somewhat 
trivial) demonstration of a non-equilibrium final state is 
a sudden quench to g = 0, where the excitation fraction 
remains constant at all times following the quench. 
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IV. TWO-PARTICLE CORRELATIONS 
A. Diffusive to ballistic crossover 

(2") 

We now consider the evolution of g s following a sud- 
den quench from g = to g = g f > 0. At t = 0, = 1 

and Vk = 0, and neglecting the quartic term g s 2 \ Eq. [4] 

simplifies to gf\t) = n 2 - gn 2 £ k e^ 5 sin ^f t} e k . At 

long times, the phase of the sin(E^t) 2 oscillates rapidly 
and we can simply replace s'm(Ekt) 2 — > (sin(-Efct) 2 ) = 



1/2. We find that g\ >(t -> oo 

,(2) 



-2(5 
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8 = 8/C,. At long times, 17^ appears to diverge as 6 — > 
as 1/5. On distances much shorter than £, one must also 
include the contribution from g$, which scales as l/S 2 . 

Therefore, at small distances, gf^ ~ l/S 2 — a/S [l9j . 

Focussing on the longer range physics, we introduce 
dimensionless variables 8 = 8/C and t = t/rMF and write 



(2) 

9s 



TT Z 8 



k 2 + 2 



(8) 



P(fc 2 + 2)1 



where k = v2k/C, and we have neglected g( 2 \ 

In Fig. [2] we plot the evolution of the density-density 
correlation function g^ — n\ as a function of t for (5 = 4, 
normalized to its asymptotic value. At short times, the 
correlations rapidly oscillate. The temporal period of 
oscillations increases with time. The structure is re- 
vealed in a saddle point approximation (with the domi- 
nant wave- vectors near k ~ 8 /t) . Within this approxima- 
tion, g^p oscillates as sin(<5 2 /t + <f>). At very long times, 
only particles with k ss contribute to the dynamics and 
the correlations saturate. As in the case of the one-body 
correlations, the requirement for correlations to saturate 
is (t — 8/c)~> tmf- 

The temporal location of the last minimum in the 
correlation function (t m i n ) is denoted by an arrow in 
Fig. [2jlcft) . On the right, we show how this minimum 
disperses with 8. When 8 < ( = h/ ^/rngfUQ, the cor- 
relations spread diffusively, while for 8 > £, correlations 
spread ballistically. This crossover arises from the nature 
of the Bogoliubov dispersion. At short distances, the 
large momentum structure of the dispersion dominates 
Ek ~ fc 2 , leading to diffusive dynamics; longer range cor- 
relations are governed by the linear part of the dispersion 
and spread ballistically. 

To further illustrate the existence of this crossover, 
we also consider a quench in the reverse direction from 
We find that gf\t) = nW(t = 



gi > to gf = 0. 



0) - 2n Q J2k e sin(e fe i) 2 M fc (0H(0). Note that even for 
a quench to the non-interacting state, density-density 
correlations evolve in time (though the condensate frac- 
tion remains fixed). As shown in Fig. fright), both 
short and long range correlations spread diffusively (we 
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where FIG. 2: Dynamics of density-density correlations in 3D 



Left: density-density correlations g ( (5/C, = 4) — n normal- 
ized to the asymptotic value at long times for a quench from 
a non-interacting state to some final gf > 0. Length and 
time is measured in terms of the condensate healing length 
(£ = h/ \jmg /no) and mean-field time Tmf = h/ (g/^o) in the 
final state. Correlations develop in an oscillatory manner and 
rapidly saturate at times t — S/c > tmf- Inset shows the 
temporal location of the minimum in the correlation function 
(arrow on left graph) (abscissa) plotted versus 8/C, (ordinate). 
Long range correlations 5/C ^> 1 spread ballistically while 
short range correlations 8/C -C 1 spread diffusively. The lines 
show purely linear and linear plus square root fit to the data. 
The latter captures the data accurately. Right: Dots are the 
same as the inset on left. Solid line shows the location of the 
minima in the correlation function for a quench from finite 
interactions to g = 0. The dynamics are purely diffusive in 
this case, with a diffusion constant ~ n. 



find 8 = TT~\/t m i n ) as the non- interacting dispersion is 
E k = k 2 /2m for all k. 

Just as with the single particle correlations, for a sud- 
den quench from g^ to gf, density-density correlations 
typically do not attain their equilibrium values at long 
times zero temperature values. For example, for a quench 
from gi to gf = 0, we find that while at long distances 

Sa>c(* = 00 ) n o ^ n n ex whereas for short distances 

^5<c(* = 00 ) ^ n o + ^non ex . Interestingly the density- 
density correlations at long times following a quench to 
a non- interacting final state do not exhibit the 1/8 diver- 
gence at small distances. In Fig. [3] we plot the density- 
density correlations at t — > 00 after a quench to gf = 0. 



B. Lattice vs. continuum 

We now ask whether additional features are found in 
the presence of an underlying lattice. The conceptually 
simplest case to consider is a quench from finite inter- 
actions to g = 0, in the presence and absence of a lat- 
tice. Here we consider a ID lattice with band structure 
e(k) = 4J sin 2 (kd/ '2), but the physics described here is 
qualitatively similar in higher dimensions. 

From Eqs.[2]and[4j we find that 
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FIG. 3: Long-time density-density correlation function 

Density-density correlation function at long times following a 
quench from finite interaction strength to g = 0. The equi- 
librium value of g 2 (5) for a non-interacting condensate is Uq. 

riQ+2non ex whereas 



Instead at long distances <?^ ? (i = oo) 



for short distances (t = oo) — > Uq +4non ex . Correlations 
fall off exponentially on a length scale proportional to the 
coherence length of the initial interacting gas. 



(2), 



(9) 



2nd 



dke 



ikS 



sin(4Jsin 2 (fc/2)t) 2 



sin(A;/2)^4Jsin 2 (fc/2) + 2 gi n 



where 5 = 5/d. Normalizing the energy and time units by 
J, we numerically integrate Eq. [S] and plot in Fig. 0] the 
typical post-quench dynamics of the density-density cor- 
relations as a function of time, following such a quench. 

At long times, the correlations in the lattice and 
continuum both decay to their asymptotic values alge- 
braically (with different exponents), but in addition, the 
lattice introduces periodic oscillations. The period of 
these oscillations is largely independent of 5, and is pro- 
portional to the time it takes for particles to travel a 
lattice unit. These oscillations are a purely lattice ef- 
fect and were also observed experimentally for quenches 
within the Mott insulating phase |(|. 

The differences between the lattice and continuum are 
even more striking at short times, as these reveal the high 
momentum structure of the underlying excitation spec- 
trum. In Fig. @] we plot the density-density correlation 
function for different values of 5 and highlight two fea- 
tures. For comparison we plot the density-density cor- 
relations in the continuum for the same quench. The 
open circles denote the position of the first maximum in 
the correlation function while the filled circles denote the 
minimum in the correlation function. On the right we 
show that these features disperse differently, the position 
of the maximum disperses ballistically, while the position 
of the minimum disperses diffusively. 

The linearly dispersing maximum has no analog in the 
continuum and emerges purely due to the lattice band- 
structure, which imposes a maximum velocity on the 
growth of correlations. In ID, this velocity is AJd [2o| . 

For quenches to an interacting final state gf ^ 0, the 
dynamics become more complicated. The minimum in 
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FIG. 4: Dynamics of density-density correlations in a 
ID lattice (Left): Typical structure of density-density cor- 
relations g^ 2 -* (<5) (t) — (/ 2 '(5)(0) for different values of 8 nor- 
malized to the asymptotic value at long times for a quench 
from gi > to gf — 0. Lengths and times are measured in 
terms of the lattice spacing and inverse hopping J -1 . Cor- 
relations in a lattice oscillate with a period proportional to 
1/J and decay with a time constant independent of 5. The 
temporal location of the first maximum (t max ) and the min- 
imum (tmin) are indicated by open and filled circles. (Right, 
Top): A plot of the evolution of the density-density correla- 
tions in the 3D continuum for the same quench at S = 5£ for 
comparison. The minimum in the correlation function (filled 
circle) disperses diffusively (Right, Bottom): t max disperses 
ballistically at all distances, while t m i n disperses diffusively. 
This latter features is the analog of the quantity highlighted 
on the top right. 



the correlation functions disperses diffusively, then bal- 
listically, similar to the 3D continuum case. The maxi- 
mum in the correlation function disperses linearly with a 
velocity set by the band structure for small values of <?/ 
crossing over to the superfluid velocity for larger gf. 



V. SHORT DISTANCE STRUCTURE OF 
TWO-PARTICLE CORRELATIONS 

Thus far we have neglected the contribution to the cor- 
relation functions arising purely from the non-condensed 
atoms. For S ~ £, and weak interactions, these terms are 
of order 0(a/() 2 , (for typical densities a ~ 50nm and 
£ ~ /mi) and can be ignored. However for 5 <C £> these 
terms give rise to a 1/S 2 divergence in the correlation 
function as a result of the singular nature of the rela- 
tive wave-function of two particles interacting with a zero 
range interaction [l9| . The two-particle correlation func- 
tion g( 2 '(5) ~ (1 — a/8) 2 where the correlations between 
non-condensed particles dominate the short distance be- 
havior and the interaction between the condensed and 
non-condensed particles leads to a sub leading 1/6 cor- 
rection which is opposite in sign. 

Recently Shina Tan showed that for a two-component 
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FIG. 5: Dynamics of the contact C(S) in a 3D Bose 
gas: (Color Online) (Left) Dynamics of the position depen- 
dent contact C(<5) = 8 2 g ( - 2 \&) (normalized to the true contact 
within the Bogoliubov approximation) at different times fol- 
lowing the quench, as a function of 8. At t < 0, the gas is 
non-interacting and C(S) — 0. Immediately after the quench, 
the zero distance correlations respond instantaneously. Col- 
ors correspond to different times: t/TMF = 0.025 (green), 0.1 
(red), 0.3 (blue), 0.9 (black). We choose a = £. For times 
t > Tmf, C(S) saturates to it's equilibrium value for a given 
8. (Right) Dependence of C(8) with scattering length a mea- 
sured in units of f. From top to bottom: 5/C = 1, 5/C = 0.1 
and 5/C = 0.01. In order to obtain the expected quadratic 
dependence of the true contact C with a for a Bose-Einstein 
condensate, one has to probe correlations on length scales 5 ~ 
nm. 

Fermi gas interacting with a contact interaction, the 
short distance structure of the two-body correlation is 
related to the internal energy via the "contact" defined 
as C = limj-vo 5 2 g^ 2 \5) [20] . For a Bose-Einstein con- 
densate, C to leading order is 167r 2 riQa 2 (l + 0(V na 3 )) 
[l9l [231 ] . Experiments using Bragg spectroscopy to probe 
the Fourier transform of g^ have a typical momentum 
resolution of k = 2ir /8 min where 6 m i n ~ 0.5/mi [2l| . As a 
result, experiments actually probe C(S) = 5 2 g( 2 '(8) eval- 
uated at Smin. 

The dominant contribution to these correlations from 

~(2) 

the non-condensed atoms is a term of the form g\ = 
|£ p e ip '*ti£(t)v p (i)| 2 . Adding this term to Eq. gj and 
multiplying by S 2 we numerically integrate the resulting 
equation to obtain C(S) for a quench from gi = to 
.9/ = 9- 

In Fig. [5]we plot C(S) as a function of 8 at different 
times. The contribution to g^ from correlations between 
the condensate and the non-condensate atoms dominates 
unless 5/C -C 1 or a > £. We choose a = 5C- Immediately 
following the quench, the zero range correlations jump to 
their equilibrium value. However C(S) changes in time for 
finite 5 and relaxes to a stationary value (denoted by the 
black curve) on times t ~ tmf. Therefore, although the 
true contact does not have any time-dependence, the ex- 
perimentally relevant quantity C(S) changes in time fol- 
lowing an interaction ramp. 

Finally we remark that even in equilibrium the scat- 
tering length dependence of C(S) is strongly influenced 
by 6. For example, the S ~ [im scale studied in experi- 



ments, one has to go to very large interactions in order 
to suppress the contribution to g^ arising from correla- 
tions between the condensed and non-condensed atoms. 
We illustrate this in Fig. [SJright) where we plot C(S, a) 
as a function of a for different values of 8, after the sys- 
tem has equilibrated following the ramp. We find that 
for 5 ~ C (typically /im), C(a) varies linearly with a at 
small a. The quadratic dependence of C(S) with a for a 
weakly interacting BEC emerges only for small 5 ~nm. 



VI. CONCLUSIONS 

Relating the information contained in the fluctuations 
of an interacting system driven out of equilibrium to the 
underlying many body parameters is a challenging task. 
Here we study a simple model, a homogenous Bose gas 
at zero temperature, where one can calculate the one and 
two body correlations of the system following a sudden 
change in the interactions. Our study reveals that the 
dynamics of these correlations reveal a wealth of infor- 
mation about the underlying excitation spectrum of the 
system. 

We first consider the dynamics of single particle corre- 
lations, such as the condensate fraction or equivalcntly, 
the excited fraction of atoms following a sudden inter- 
action ramp. For a quench from a non-interacting gas, 
we find that all the dynamics are governed by a single 
timescale, the mean-field time. For quenches between 
two interacting systems, we find two timescales governing 
the long and short term dynamics of the gas. Moreover 
we find that the excited fraction at long times after the 
quench is always greater than the expected equilibrium 
value. 

We then consider the dynamics of the two particle cor- 
relation function. For quenches between interacting ini- 
tial and final states, we find a non-trivial crossover be- 
tween diffusive spreading of short-range correlations and 
the ballistic spread of long range correlations. We relate 
this crossover to the underlying Bogoliubov dispersion. 
We then consider the case of a Bose gas in a lattice and 
discuss the additional features that arise in the correla- 
tions from the band structure. Notably, we find that the 
band structure leads to an additional linearly dispersing 
feature in the correlation functions. 

Finally we discuss the dynamics of the contact follow- 
ing a sudden quench. We find that while the true zero 
range contact instantaneously takes on the new equilib- 
rium value, the finite resolution of an experiment will 
make the contact appear time dependent. 
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